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In this manuscript we present a comprehensive study on the multifractal properties of high- 
frequency price fluctuations and instantaneous volatility of the equities that compose Dow Jones 
Industrial Average. The analysis consists about quantification of dependence and non-Gaussianity 
on the multifractal character of financial quantities. Our results point out an equivalent influence 
of dependence and non-Gaussianity on the multifractality of time series. Moreover, we analyse i- 
diagrams of price fluctuations. In the latter case, we show that the fractal dimension of these maps 
is basically independent of the lag between price fluctuations that we assume. 
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I. INTRODUCTION 



Scale invariance and fractality, i.e., the absence of a characteristic scale can be found in a widespread of natural 
and human creation phenomena Mathematically, scale invariance of a certain function, /, of an observale O is 
written as, 

/ (XO) = A" / (O) , (1) 

and it has consensually been considered as a signature of complexity Q • Specifically, power-law behaviour respecting 
Eq. ([I]) has empirically been verified in the probability density and autocorrelation functions of several time series such 
as: fluctuation in heart rate beating [3 .gait [|, variation in the magnetic field of the solar wind in the hclioshcath 
Q, or relative stock price fluctuations Q0] among many others. Concerning time series and fractality 1 , if many 
of them seem to be monofractal Q, i.e., they are characterised by a single scale exponent, just as in Eq.dU), other 
time series, namely those we have referred to here above, have shown a spectrum of locally dependent a exponents. 
Analytically, this is noted as 

f({XO} v ) = X a(v) f(O v ). (2) 

The previous Eq. |(5J) is also valid for multiscaling and multifractality as well, which has consistently been associated 
with the main statistical features of time series obtained from complex systems. Consequently, this close relation has 
been prominent in either the development of dynamical models or validation of previous approaches. In the former 
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case, pioneering works by B. Mandelbrot have opened the door to a new treatment of financial markets dynamics 
9]. 

In sequel of this manuscript we perform an extensive analysis of the statistical features of high-frequency price 
fluctuations, r t , of the 30 equities that compose the Dow Jones Industrial Average (DJIA). Previous studies on daily 
price fluctuations have shown the existence of a multifractal behaviour [icl |. Hence, with this high-frequency analysis, 
it is our aim to study the multiscaling of price fluctuations at a level that is closer to the transaction dynamics as 
it has been made for other financial observables. Our study is driven on the evaluation of the multifractal spectra 
of both of time series and (£ = l)-diagrams, (r t ,r t+ i), describing their main factors of multifractality. In addition, 
we enquire into price fluctuations absolute values, |r t |, also called as instantaneous volatility, vt = \r t \, multifractal 
behaviour and analyse its weight on the multiscaling characteristics of price fluctuations. Our manuscript is organised 
as follows: in Sec. [Til we describe the data used and the methodology applied in order to obtained multifractal spectra 
for time series and ^-diagrams. Along Sec. [TJTJwe present our results of the analysis of price fluctuations and volatility 
multifractal spectra. This comprises the quantification of the key elements of multiscaling for both quantities. In 
addition, we verify the plausibility of a superstatistical approach (which is able to provide a nice answer within the 
context of price fluctuations probability density function) in a multifractal characterisation of price fluctuations. In 
Sec. II VI we present the results of the study of the fractal dimension of price fluctuations ^-diagrams. To finalise, some 
remarks, conclusions, and perspectives for future work are set forth in Sec. [V] 



II. DATA AND METHODS 



A. Data 



Our data is composed by 1 minute time series of the prices, Si (t) (i stands for the company), of the 30 companies 
that composed the Dow Jones Industrial Average from the 1 st of July until the 31 s * December of 2004 in a total of 
around 50 000 data points for each equity. For each i equity we have computed 1 minute (log-)price fluctuations as, 

f'i (t) = In Si (t)- In Si (t-1) (3) 

It is well known that trading activity exhibits a intraday pattern [11] . In other words, markets tend to be highly active 
(hence volatile) in the first 30 minutes of a business day, mainly to take advantage from news and events between 
the closure of the previous market session and the next following opening. After a decrease of activity along the day, 
markets present an activity set-up in the final part of trading sessions, basically due to the action of liquid traders. 
This U-shape enhances spurious features namely in correlations. To remove it, we have performed according to the 
following standard procedure: 

• After we have computed 1 minute price fluctuations, as in Eq. @, we have determined the average volatilities, 
A, associated with equity i and intra-day time t' (which has an upper bound of 340 minutes for companies 
traded at NYSE and 420 minutes for companies traded at NASDAQ), 

(4) 

where N represents the number of days for which the market was trading at t' intra-day time; 

• We have then used the average volatilities to normalise price fluctuations, eliminating the intraday pattern, 

h (f) , \ 

^aTW < 5 > 

where we have dropped the prime of t in the left-hand side, because the time series has lost its intra-day profile; 

• To complete, we have removed the average and normalised {fj (t)} by its standard deviation, 

* (*) = fiMzM , (6) 

((h) 2 ) -(h) 2 

((. . .) represents average over all the elements of time series i), to define our studied price fluctuations, Ti (t). 



We next present our methods to evaluate multifractal spectra for time series, multifractal detrended fluctuation 
analysis (MF-DFA), and ^-diagrams. 



B. MF-DFA 



The MF-DFA I12j is one of the most applied methods to determine the multifractal properties of time series in 
several fields [H| [lTj • We have chosen to apply MF-DFA in lieu of Wavelet Transform Modulus Maxima (WTMM) 
[l8[ taking into account a recent comparative study where it has been shown that in the majority of situations MF- 
DFA presents reliable results [191 ], i.e., it does not introduce specious multifractality, at least in the amounts that have 
been computed from WTMM. The MF-DFA method goes as follows: 

Consider the time series {xi (t)} (xi represents both of price fluctuations and instantaneous volatilities of company 
i) composed by N {N > 1), 

• Determine the profile Yi (t) that corresponds to the deviation of signal elements from the mean 

t 

Yi (t) = Y, fa (*) - <^>] ,(!<*<#), (7) 
i=i 

and thereafter, 

t 

Y t (t) = Y / [Y t (t)-(Y i }}; (8) 
l=i 



• Divide profile Yi (t) into N s = int (— ) non-overlapping intervals of equal size s; 

• Compute local tendency by a least-square adjustment, and thereupon variance 2 , 



F 2 (M = [{y-l)s + l]-y v {l)\ 

i=i 



(9) 



for each segment u, v — I,.. , ,N S , where y v (i) represents a m'^-order polynomial. The order of the polynomial 
is relevant on the results one might obtain. For the series we have analysed we have used polynomials of order 
5 from which on we could not appraise changes of the values of multifractal spectra. 



• Figure out the average F z (s) over all segments to obtain the fluctuation function of order z, 

V* 



r i J^r. wO 



V 



and 



F z (s) = exp 



z = 0. 



(10) 



(11) 



• Assess the scaling behaviour of F z (s) considering log — log scale representation of F z (s) vs. s for each value of 
z. In case the series {x (i)} shows multiscaling features then, 

^i£)~ s M^). (12) 



Small fluctuations are generally characterised by large scale values of exponent h(z) (and z < 0), whereas large 
fluctuations are typified by small values of h (z) (and z > 0). 

To bridge this procedure with the standard formalism, we can verify that [F z (s)] can be interpreted as the partition 
function, Z z (s) 8], which is known to scale with the size of the interval as, 

Z z (s)^s T ^. (13) 



2 In the rest of this section we omit company index i to turn out notation lighter. 



Hence, according to Eq. (|12p and Eq. (|T5|) we have, 

r(z)=zh(z)-l. (14) 

Using Legendre transform, 

f(a) = za-r(z), (15) 

we can relate exponent r (z) with Holder exponent, a, 

, . . dh (z) , . 

a = ft (16 

and 

/ (a) = z [a -ft (?)] + !• (17) 

For z = 2, ft (2) = _ff , that corresponds to the Hurst exponent [20l | customarily determined by methods like the original 
R/S ratio or DFA [2l| (for nonstationary signals) from which MF-DFA derives. For z = 0, / (a) obtained from Eq. 
(fT5")) . Eq. (fl7o|) , and Eq. (fl"T|) corresponds to the support dimension. 

In the case of a monofractal, h (z) is independent from z, since there is homogeneity in the scaling behaviour. 
Specifically, there exist only different values of h (z), for each z, if large and small fluctuations scale in different ways. 

C. Box counting algorithm 

Box-counting methods have been extensively applied to determine scaling (fractal) properties of measures. As a 
matter of fact, it is the standard procedure to verify the factal nature of measures p. To determine multiscaling 
properties in the ^-diagrams we have used a recently presented optimised implementation [15j of the procedure 
introduced in Ref . [221 ] . This method is an improvement of the algorithm proposed by Liebovitch and Toth [23j , and 
it is based on the fact that coordinates of a fractal, suitably shifted and rescaled, and written on the binary numerical 
base can be combined to form bit strings with k De bits whose first itiDe bits from left to right determine uniquely 
the position of the coordinates in D^-dimensional space. Here, m = 1, 2, • • ■ , k and k is a positive integer. Thus, the 
N De coordinates are mapped in N bit strings and with kDE bits, where k is the maximal number of bits used to 
represent each coordinate on the binary base. After masking uiDe bits from right to left, strings that have the same 
position code belong to the same box in resolution m. Then, by scanning the N bit strings, the number of changes are 
stored. The number of changes represents the number of boxes needed to cover the fractal set in the scale m = ln2 s. 
If the set presents a fractal measure, one has 

N~s~ D f, (18) 

where Df is the fractal dimension of the set. 



D. Quantification of multiscaling components 

As we have stated in Sec. Q] there has been established a close relation between multiscaling and both of correla- 
tions and probability density functions of price fluctuations time series. In another perspective, we can ascribe to 
memory and non-Gaussianity of probability density functions the emergence of multifractal characteristics in price 
fluctuations This can be made if we first consider such contributions as independent. Upon this assumption, 

we can quantify their relative weights. In other words, if we aim to size up the weight of non-Gaussianity we must 
destroy memory in the signal. And from it, by using the independence conjecture, we determine memory influence. 
On one hand, memory is basically destroyed if we shuffle time scries elements. Doing that, we reorder the values of 
our original time series, but we keep the stationary probability density function unalterable. On the other hand, we 
can destroy non-Gaussianity by implementing the procedure which we call as phase randomisation: 

• Determine the Fourier transform of the signal {x {t)}, 



£ f = F[x(t)]; 



(19) 



• Dissociate amplitude from the phase of the transformed signal, 



£/ = 10 1 ex P 



i arctan 



where Im(Re) stands for imaginary (real) part of some complex number 



(20) 



• Expunge the phase of the transformed signal and introduce new random phases, 9 rn d, uniformly distributed, 
for half of the elements, and assign for the other half of the series a phase —Q rn d- In this way we define a phase 
randomised signal, 



|£/|exp[z 9 rnd ] ; 

• Apply the inverse Fourier transform on the phase randomised signal, 



x rnd (i) = JF" 1 [C f nd ] 



(21) 



(22) 



For both shuffled and phase randomised time series obtained from the original signal, we can also carry out a 
MF-DFA analysis. For each case, Eq. (|12j) can be verified where we use exponents h s hf ( z ) for the shuffled time series, 
and h rn d (z) for the phase randomised case. Assuming independency between multifractal factors, we have measured 
the contribution of correlations, h cor (z), by, 

h cor (z) = h (z) - h s hf (z) ■ (23) 

If only these two factors introduce multiscaling on the signal then, when we perform the phase randomisation process 
on a shuffled signal, we should obtain a Gaussian and uncorrelated signal, z.e., h s }if— rn& 

(z) = i for all z. Theoretically, 
we can evaluate the contribution of non-Gaussianity, h' PDF (z), from phase randomised time series as well, 

h' PDF (z) = h sh f (z) - hshf-md (z) ■ (24) 

However, the probability density function of a finite time series is influenced by its size, particularly for small time 
series In this sense, comparing results obtained from times series with different probability density functions, 

such is the case of {x sh f (i)} and { x sh f~ rnd (t)\, introduces error factors that we are not able to quantify. Regarding 
this factor, we have opted to define an effective contribution of non-Gaussianity, hpuF (z), 

hpDF (z) = hshf (z) ■ (25) 

In a previous article by us [16| |. in which we analyse the multifractal features of traded volume for the same equities, 
we have computed a contribution that we have called as non-linear effects. These effects are actually governed by 
finite size effects that play a significant role on the multifractal character of a time series when large values of \z\ 
are taken into account. Moreover, the finiteness of a time series might introduce fake multiscaling features. This 
fact emphasises the sensitiveness of multifractal measurements which are many time inflated by artefacts. In order 
to avoid, or at least minimise those spurious features, a careful choice of the range of s and z values must be made. 
In our analysis, we have chosen s between 8 and 11585, and z between —3 and 5. Within this range of values we 
were able to obtain numerical curves which concur to the theoretical scaling curve of independent and Gaussian time 
series. 

Multifractality can be effectually quantified through the difference between scale exponents of z m ; n and z mSLX , 

Ah = h (z min ) - h (z max ) . (26) 

For a monofractal, Ah = 0, because of the linear dependence of r with z. Equation (p?6"|) can be used for the 
original time series, Ah, and for the shuffled time series, Ah s hf. From these values, we finally compute the weight of 
non-Gaussianity, Ah s hf / Ah, and of correlations 1 — Ah s hf / Ah. 



III. RESULTS FOR TIME SERIES 



A. Multifractality for price fluctuations time series 



In Fig. Q] we present our results for the multifractal spectrum of price fluctuations, shuffled, phase randomised, 
and shuffled plus phase randomised time series. The values have been obtained by performing an average over the 



30 companies of moments r. Despite of the fact that it has been verified the influence of equities liquidity on the 
multifractal properties of financial time series, all companies of our data set have presented liquidity values within the 
same order of magnitude, turning out our average over the companies perfectly plausible. As it can be seen, the price 
fluctuations time series present a wide multifractal spectrum with a m i n = 0.364 and a max = 0.724. Furthermore, we 
verify a strong asymmetry between the part of the spectrum that goes from a m j n up to a (z = 0) and the remaining 
part of the spectrum. The asymmetry in Fig. [T] is contrary to the / (a) curve that has been measured in fully 
developed turbulent flows [24| often considered a price fluctuation analogue. Concerning the other time series, we 
observe that the multifractal spectrum of the shuffled time series is less wider than the spectrum for the original time 
series. In addition, the shuffled signals have larger spectrum than the randomised and shuffled plus phase randomised 
signals. Analysing scaling exponent h(2), that is the common Hurst exponent, H, we have obtained a value around 
\ concomitant with a white noise sequence, and in accordance with the Efficient Market Hypothesis (EMH) [25| . 
Furthermore, we have max (/(a)) = 1, i.e., price fluctuations time series are fat-fractals as it occurs for a large 
variety of other signals and non-linear phenomena [l[ [H[ . For the h difference defined in Eq. (|2"6"|) we have obtained 
Ah = 0.15 and Ah s hf — 0.08. These values yield a weight of 54% for non-Gaussianity and 46% for correlations in 
the multifractal properties of our time series. In spite of this result appears to be at odds with the H = 4, we must 
call attention to the fact that there is a more delicate relation for random variables, the statistical dependence [27| . 
which cannot be described by the Hurst exponent. The statistical dependence of financial observables [(| (28[ [2!| has 
been verified by means of mutual information measures (30l | . We attribute to this statistical feature the multiscaling 
of price fluctuations we have perceived. This assignment is also supported by the structure of {I = l)-diagrams that 
we analyse in Sec. IIVI 




FIG. 1: Multifractal spectra / vs. a of the price fluctuations (■), shuffled time series (•), phase randomised (♦), and shuffled 
plus phase randomised (A) time series of DJIA equities. As it is depicted, when elements that introduce multiscaling are 
removed the multifractal spectrum becomes narrower. 

In Fig. [5] we show over different panels the moment r as a function of z. We observe that only the shuffled plus 
phase randomised signal are in compliance with the theoretical curve, r = z/2 — 1, of an Gaussian time series of 
independent elements. 



FIG. 2: Scaling exponent r vs. z of the price fluctuations (upper left), shuffled (upper right), phase randomised (lower left), 
and shuffled plus phase randomised (lower right) time series of DJIA equities. In all panels, the line r = | — 1 represents the 
theoretical curve for an independent and Gaussian time series. For the shuffled plus phase randomised time series, the points 
coincide with the curve of a Gaussian time series of independent elements. For the other cases, we verify a departure from the 
line, proving the multifractal behaviour of price fluctuations. 

B. Multifractality for instantaneous volatility time series 

Albeit volatility is not directly observable, it plays a central role in financial modelling [3l|, and it is usually related 
to the magnitude of price fluctuations. It is on this quantity that long-lasting covariances associated with asymptotic 
power-laws are measured. As a matter of fact, the appropriate mimicry of a long-lasting autocorrelation function of 
the volatility associated with a white noise character of the variable upon study is one of prime challenges in several 
areas of scientific research. Aiming to appraise its potential multiscaling nature we have performed a MF-DFA analysis 
on instantaneous volatility time series. The main results are shown in Fig. [3] and Fig. 2] From our analysis, we have 
verified that there are clear differences between multifractal spectra for price fluctuations and absolute values. 

In first place, and against our primary expectations, we have observed that price fluctuations have a wide multifrac- 
tal spectrum. Specifically, we have computed Ah — 0.15 for price fluctuations, and Ah = 0.10 for volatilities. This 
corresponds to a ratio of 2 over 3. As it happens for price fluctuations, the multifractal spectrum is asymmetric. We 
have also obtained h (2) = 0.71, which indicates a strong persistency on volatility time series in accordance with previ- 
ous empirical findings. We clarify that we expected to obtain a wider spectrum for instantaneous volatility because of 
correlations and non-Gaussianity of this quantity. For shuffled instantaneous volatility time series we observe a shift of 
/ (a), and a lessen of curve width. On the other hand, when we turn instantaneous volatility into a Gaussian variable 
points multifractal tends do be clearly diminished, though still present. This is in accordance to previous verifications 
about local fluctuations on Hurst exponent for financial time series [45| which introduce multifractality. Bearing in 
mind the value Ah = 0.10, the difference between scaling exponents of the shuffled time series, Ah s hf — 0.05, points 
non-Gaussianity and dependence as equally responsible for the multiscaling of instantaneous volatility. From Fig. 2J 
it is visible that r s hf almost coincides with the theoretical curve of an independent and Gaussian time series. Such a 
result indicates that the probability density function presents a nearly exponential decay. We corroborate this result 
with Fig. H] in which we present absolute values probability density function, p(v). In line with Fig. [5] we verify that 




FIG. 3: Multifractal spectra / vs. a of the instantaneous volatility (■), shuffled time series (•), phase randomised (♦), and 
shuffled plus phase randomised (a) time series of DJIA equities. As it is shown, when elements that introduce multi-scaling 
are removed the multi-fractal spectrum becomes clearly narrower. 



p (v) fits for a F-distribution, 



/"u\'l > r v 



(27) 



where 9 = 0.32 ± 0.02, <fi = 1.83 ± 0.01, and q — 1.08 ± 0.02. Taking into account error margins, the small deviation 
from exponential decay given by numerical adjustment is in agreement with the slight deviation of r s hf from the 
theoretical curve that we have measured. 



C. Effects of the signal and (instantaneous) volatility multifractal behaviour on price Actuations multiscaling 

In this subsection we assess the influence of the multifractal character of instantaneous volatility on the multifractal 
nature of price fluctuations. To that, we have proceeded the following way. We have separated price fluctuations, 
r (t), considering each element as the product of elements of two other time series, i.e., one that considers the signal 
of the price fluctuation, s (t) = ±1, and other which takes into account the magnitude or instantaneous volatility, 
v (t) — \r(t)\. Preserving the signal time series, we have multiplied {s (t)} by time series that were obtained after 
shuffle, v s hf{t) 1 phase randomisation, v rn d(t), and shuffle plus phase randomisation, v s hf-md{t), procedures. The 
results we have obtained are depicted in Fig. El and Fig. [7] From Fig. El we see that the statistical properties of 
volatility do influence the multifractal spectrum of price fluctuations. If we only shuffle {v (t)} elements, the u (t) 
time series, 

u (t) = s (t) v (t) , 

just has a paltry narrower f (a) curve than {r(t)}. It has Ah = 0.13 in opposition to Ah = 0.15 of {r(t)}. 
This is an unexpected result regarding the influence of {v (t)} ordering on its multifractal spectrum. However, 
when we destroy the non-Gaussianity of instantaneous volatility probability density function, we basically destroy 
the multifractal spectrum of price fluctuations, since Ah = 0.04, or Ah — 0.03 when we combine shuffling with 
phase randomisation procedures on {v(t)}. The latter result also sets the influence of the signal ordering on the 
price fluctuations multifractal character at the order of error in absolute accordance with previous analysis for other 
characteristics, namely the approach to the Gaussian when of cummulative price fluctuations probability density 
functions [4p] |. 
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FIG. 4: Scaling exponent r vs. z of the instantaneous volatility (upper left), shuffled (upper right), phase randomised (lower 
left), and shuffled plus phase randomised (lower right) time series of DJIA equities. In all panels, the solid line r = § — 1 
represents the theoretical curve of Gaussian time series of independent elements. Regarding instantaneous volatility results, it 
is visible the departure from Gaussian independent behaviour that persists when we destroy the Gaussianity. In the lower left 
panel the dotted line represents the monofractal curve r = Hz — 1 with H — h(2) = 0.71 ± 0.01. If we considered the phase 
randomised time series as a pure monofractal set we would have the best fit for H — 0.692 ± 0.002, a bit ouside error margin 
of h(2). 



As it has been observed @ [2!| , many of the dynamical and statistical properties of price fluctuations depend on the 
volatility. Although it is a pivotal variable in finance the truth is that volatility definition is still ambiguous [321 ] . If 
in many situations it is presented has we have been doing, volatility is oftenly determined as the standard deviation 
of price fluctuations over window of length I 3 . The latter definition is widely applied on stochastic volatility models. 
In that particular case, superstatistical models have been applied in problems of financial origin [6] [33] to define such 
models. Concisely, super statistics or "statistics of statistics" (3~ij is a compound method which has emerged within 
statistical mechanics. It is based on the assumption of a local statistics dependent on a parameter that fluctuates 
(smoothly) on a time scale that is very large when compared with the time needed for a system to reach a local 
equilibrium or stationarity. In a superstatistical context it has been proved that, if we have a set of local Gaussian 
random variables, 



Q a {x) 



i 



27TO- 



exp 



x 

2^2 



(28) 



and the inverse variance, a 2 , is associated with a T-distribution, 





" X' 


(f) 7exp 


~1. 



(29) 



then, the stationary distribution given by 



p(x) = J Ga(x)T(a- 2 ) d(a- 



3 When 1 = 1 we obtain the instantaneous volatility definition. 




FIG. 5: Instantaneous volatility probability density function P (v) vs. v averaged over D JIA equities. Symbols are the empirical 
PDF and the line the best fit using a F-distribution, Eq. (J2TJ) (x 2 /n = 4.4 x 10~ 6 and R 2 = 0.999). 



is equal to a Tsallis (or Student t-) distribution [35j . 



p{x) = 



where 



l-(l-«) 



(30) 



9 = 1 



27 



(31) 



In this way, superstatistics has been considered has the first dynamical foundation for non-extensive framework [361 ] 
that has non-additive entropy, S q [13], as its cornerstone. Distribution ([3"Uj) has regularly been used to fit for price 
fluctuations of several financial markets, and also for the data set we have been analysing for which it has been found 
a value of q = 1.31 ± 0.02 [l7| • If we assume a superstatistical approach for the data set upon analysis from Eq. (|3"Tj) 
we obtain 7 = 1.82. 

In what follows we analyse a discrete ARCH-like process [38| that can be catalogued as superstatistical. Explicitly, 
we have generated time series, {y (£)}, from the product of an uncorrelated Gaussian signal, {to (i)}, with (ui (t)) = 0, 

and (lj (t) 2 \ — 1 by an uncorrelated volatility signal, {a (t)}, 



y{t) = a(t) w(t), 

such that (j~ follows a T-distribution with 7 = 1.82 as we have obtained. In this case, because we neglect memory 
on volatility, we can compare the multifractal study of this time series with the results that we have presented at 
the beginning of this subsection IIII CI for {u(t)} with a shuffled instantaneous volatility. We have opted for this 
comparison because, just as s(t), u)(t) does not contribute to the multifractal spectrum. The excerpt of the time 
series we have generated is presented in Fig. [S] In the same figure we comprove that {y (t)} follows PDF f3"U]) with 
q = 1.3. 

Afterwards, we have performed a multifractal analysis along the same lines we have made for price fluctuations. 
Even though both multifractal spectra are very similar we can verify a noticeable difference. As a matter of fact 



1.0 
0.9H 



0.8- 



0.7- 



0.6 




0.3 



FIG. 6: Multifractal spectra / vs. a of the price fluctuations (■), and of time series {r (t)} that use shuffled (o), phase 
randomised (0), and shuffled plus phase randomised (A) volatility time series of DJIA equities. As it is depicted, the multifractal 
character of volatility plays an essential role at the multifractal nature of price fluctuations. This role is clear for the non- 
Gaussianity of v (t) . 




FIG. 7: Scaling exponent r vs. z of time series {r (t)} that use shuffled (right), phase randomised (centre), and shuffled plus 
phase randomised (left) volatility time series of DJIA equities. In all panels, the line r = | — 1 represents the theoretical curve 
for an independent and Gaussian time series. The importance of the multifractal characteristics of volatility are demonstrated 
by the clear approach of these three results towards the theoretical curve of an independent Gaussian process. 



we have obtained Act = 0.28 for {y (t)} with shuffled volatility time series, and Act = 0.24 for the generated series 
that we have priory analysed, i.e., an error of 17%, see in Fig. [U This means that superstatistics can be considered 
as an acceptable first approach, although models that consider long term memory in variance [39j are certainly 
more appropriate. Since the only source of multifractality in this case is the asymptotic power-law behaviour of the 
stationary PDF P (y), {y (t)} should be in fact called a bifractal with r (z) = for z > 5.45. 



IV. RESULTS FOR /-DIAGRAMS 



A rich and interesting way of representing time series is to consider a mapping of the time series onto a plane where 
each point signalled is obtained by pairing elements xt and x t +i of the time series as ordinate and abcissa, respectively. 
This ^-diagrams and related methods [43| are frequently used on studies about biological [42| , and dynamical systems 
p3 |. Moreover, they have also been introduced to study daily fluctuations of some securities [44|. This type of 




FIG. 8: Left: Excerpt of the ARCH-like signal (lower panel) where y' represents y divided by the average value of a, 
a m . The elements of a signal (upper panel) follow a PDF such that P (c -2 ) is a T-distribution with 7 = 1.82, and 5 — 2. 
Right: Stationary PDF of y' , P (y') vs. y . Symbols have been obtained from the time series and the line is the numerical 
adjustment for a (/-Gaussian distribution with q = 1.3. Although this is not an exact approach, the adjustment is rather nice 
(x 2 /n = 1.2 x 1(T 5 and R 2 = 0.99). 




FIG. 9: Left panel: Multifractal spectra / vs. a of the {y (t)} (•), and of time series {r (t)} that use shuffled (o) volatility 
time series of DJIA equities. Right panel: Scaling exponent r vs. z of {y (i)} (•) and of time series {r (t)} that use shuffled (o) 
volatility time series of DJIA equities. The line r = | — 1 represents the theoretical curve for an independent and Gaussian 
time series. 

representation, full called as l-diagram variability method [42l |. is in fact quite illustrative since it is a simple way of 
capturing regular aspects of systems which are apparently irregular. Such regularities can be characterised by regions 
which are more visited in space Xt X Xt+e- Specifically, taking into account price fluctuations time series and £ — 1 
as an example, it allows one to verify how prices evolve in segments of two time intervals. Nextly, we analyse the 
first return map of the price fluctuations. In Fig. [10] we show the plot of r t +i versus r t for some of the companies of 
our set 4 . The plots present a very interesting structure. Over the four quadrants (anticlockwise) we have got stripes 
with high density of points and "forbidden" regions close to the axes. We assign to transaction costs the emergence 
of this banned regions. In the 3 rd quadrant we can see a highly visited region close to the origin, point that small 
decreases induce small decreases. We have investigated the probabilities for each quadrant and we have found a very 
peculiar behaviour for DJ30 5 . The probabilities for each quadrant (anticlockwise) can be interpreted as follows: 



4 The plotted companies have been chosen in order to represent different sectors of activity and ways of trading (NYSE and NASDAQ). 

5 We have also calculated these probabilities for original data - intra-day trend mask the effects observed in first return maps - and, once 
again, we have found the same behaviour, but with different probabilities. 




FIG. 10: Recurrence maps (step s = 1) for the companies Caterpilar (A), Citibank (B), Intel (C) and 3M (D) - detrended 
data.(i = t) 



1 st quadrant - probability of two consecutive profits, Pi; 2 nd quadrant - probability of a profit after a loss, P2; 3 rd 
quadrant - probability of two consecutive loss, P3; 4 th quadrant - probability of a loss after profit, P4. These results 
are shown in Table HI As it is easily observable the dynamics of the system is basically up- and- down- and- up since 
the fourth and second quadrants together represent 2/3 of the points plotted on 1-diagrams, with the probabilities 
of having either two consecutive profits or two consecutive losses equal to 17%, in average. Interestingly, we have 
verified that although the number of negative price fluctuations surpasses the number of positive price fluctuations, 
N (ri > 0) — N (r, < 0) = —199 (related to the skewness of the distributions), the cumulative sum of price fluctuations 
yields a positive value for all equities, ^ t r.i (t) — 3.82 (in average). In other words, although during the period upon 
analysis there was a larger number of negative price fluctuations than positive price fluctuations, the magnitude of 
the latter were greater so that a positive evolution arose. As a matter of fact, during this period the DJIA index 
increased its magnitude from 10334.16 to 10783.01, or a heighten of 4.8%. 

Be aware that, looking at Fig. [TU1 there exists a clear pattern for these probabilities. In order to further show that 
these characteristic patterns go beyond the uncorrelated essence of price fluctuations time series, we have performed 
immediate 1-diagrams for the shuffled signals. The results are presented in Fig. [IT] where it is visible that these 
diagrams are different from the diagrams that we have shown in Fig. 1101 namely the accumulation around lines 
r t +\ = ±r t becomes less clear. Furthermore, analysing shuffled plus randomised times series, Fig. 112) we have 
observed the lost of any pattern, forbidden stripes inclusive. Actually, both of the two latter representations are more 
homogeneous. In our opinion this is a clear evidence about the importance of dependencies and non-Gaussianity on 
price fluctuations dynamics. At this point it is absolutely necessary to stress that this profile for 1-diagrams does not 
contradics the EMH, if one tried to make use of this property for immediate trading, transaction costs would surpass 
any possible (read likely) income. 

Analysing ^-diagrams for I = 2,4, 10 we have verified an equal occupancy of all quadrants, Pi = P2 = Ps = P4 = 
25%, which indicates the loss of any predictability on the time series. 

To quantify the properties of ^-diagrams we have also applied the algorithm described in Sec. Ill CI Namely, we 
have mapped the space r t x r t +\ onto interval [2°, 2 16 — 1], and we have estimated the fractal dimension of this space 
structure for the 30 companies. In Fig. [T3]it is noticeable that for the majority of companies the scale regime holds 
over a large range of scales. We have then used the interval 2 2 to 2 8 to numerically obtain the fractal dimensions that 
are shown in Table HH which correspond to the slopes of the fitting straight lines in Fig. 1131 Therein, it is verifiable 




FIG. 11: Recurrence maps (step s = 1) for the companies Caterpilar (A), Citibank (B), Intel (C) and 3M (D) - detrended and 
shuffled data.(i = t) 



that the fractal dimension vary slight as the step (s) is changed, as well as after a shuffling procedure (keeping the 
order of the diagram). However, it is strongly affected by phase randomization, and it presents for this case values 
that are compatible with a 2-dimensional Gaussian distribution. 



V. FINAL REMARKS 



To summarise, in this manuscript we have made an exhaustive analysis of the effective multifractal properties 
of high-frequency price fluctuations and instantaneous volatility of the equities that compose Dow Jones Industrial 
Average. This analysis has comprised the quantification of dependence and non-Gaussianity on the multifractal 
character of price fluctuations and volatility. Furthermore, we have studied the multifractal properties of the l- 
diagrams made from price fluctuations time series. Our results indicate that dependence and non-Gaussianity have 
similar weights on the multifractal features of both financial quantities. Contrarily to some stylised facts, and especially 
for instantaneous volatility [451 ] . we have not verified a solid asymptotic power-law decay of the probability density 
functions, i.e., fair deviations from exponential decay. This result is substantiated by the clear approach of r (z) 
curves to the theoretical curve of a independent gaussian signal when we perform a shuffling on time series elements. 
If we consider persistence as a major factor for multiscaling, it might be puzzling to verify that multifractality for 
price fluctuations is stronger than it is for magnitude price fluctuations. Such an apparent contradiction is cleared up 
if we take into consideration that price fluctuations PDF appears to be more fat tailed than instantaneous volatility 
which introduces a larger contribution to multiscaling. Besides, in respect of probability density functions, we have 
observed that a superstatistical approach to price fluctuations appears to be valid as a first approach. Still on 
multiscaling, we have tried to appraise the robustness of instantaneous volatility by means of measuring the effect of 
its possible multifractal nature on price fluctuations multifractal properties. Our results have indicated that the non- 
Gaussinity of instantaneous volatility (price fluctuation magnitudes) is the chief element of multifractal properties of 
price fluctuations. This occurs because the uncorrelated character of the signal annihilates the influence of dependences 
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FIG. 12: Recurrence maps (step s = 1) for the companies Caterpilar (A), Citibank (B), Intel (C) and 3M (D) - detrended, 
shuffled plus phase randomised data.(i = i) 



of instantaneous volatility leading to the non-Gaussianity of latter quantity the chief role of introducing multifractality 
on price fluctuations time series. In this perspective heteroskedastic (i.e., ARCH) approaches, within superstatistics 
is enclosed, to price fluctuations are validated. 

Analysing ^-diagrams obtained from price fluctuations time series we have got sequences of immediate price fluc- 
tuations around Cartesian axes that are forbidden. We have attributed this fact to transaction costs. We have 
also observed that despite the number of negative price fluctuations is greater than the number of positive price 
fluctuations, the sum all returns is in fact positive, which is in accordance with both price fluctuations skewncss 
and economical evolution. By means of a box counting algorithm we have computed the fractal dimension of such 
diagrams. We have verified that the fractal dimension varies slightly when time ordering is destroyed, and it is deeply 
affected by randomisation procedures. This provides an important clue on the fundamental role of non-Gaussinity of 
price fluctuations in several properties usually observed. 
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TABLE I: Probabilities for each quadrant (columns 2-4) of 30 companies of the DJ30. In column 6 is shown the difference 
between positives price fluctuations and negative price fluctuations, and in column 7 is shown the sum of all returns. The last 
two columns have been obtained from trended data. Even though most price fluctuations are negative (for most companies) 
though the sum is positive. Note also that there is a clear pattern for the quadrants for the 30 companies. 
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